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Abstract 

A certain notion of canonical equivalence in quantum me- 
chanics is proposed. It is used to relate quantal systems 
with discrete ones. Discrete systems canonically equivalent 
to the celebrated harmonic oscillator as well as the quar- 
tic and the quasi-exact ly- solvable anharmonic oscillators are 
found. They can be viewed as a translation-covariant dis- 
cretization of the (an)harmonic oscillator preserving isospec- 
trality. The notion of the g— deformation of the canonical 
equivalence leading to a dilatation-covariant discretization 
preserving polynomiality of eigenfunctions is also presented. 
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Two classical mechanical systems related by a canonical transforma- 
tion are equivalent (for general discussion see, for example, [Q). It is 
quite natural to try to adopt a similar notion for the quantum mechani- 
cal systems making a difference that classical canonical transformations 
should be replaced by their quantum counterpart, which means that 
the Poisson bracket is replaced to the Lie bracket, [P{p, q), Q{p, q)] = 1- 
Present work is an attempt to explore this natural definition taking 
as an example both (an) harmonic oscillator and a certain particular 
type of quantum canonical transformations. There are well-known dif- 
ficulties on this way and we will try to indicate and overtake them. 
Meanwhile, we will call systems related through a quantum canonical 
transformation canonically equivalent. 

Let us remind that any quantum-mechanical system with a Hamil- 
tonian 7i{x,p) is intrinsically related to a Heisenberg algebra [x,p] = 
—i. Hamiltonian can be treated as an operator (or, an element of 
the Heisenberg- Weyl algebra) acting in the quantum-mechanical phase 
space which is defined as an object consisting of the universal envelop- 
ing Heisenberg algebra (the Heisenberg- Weyl algebra) and the vacuum 
|0 >, such that p|0 >= 0. From this point of view, canonical transfor- 
mations are nothing but changes of variables in the phase space pre- 
serving polarization or, in other words, preserving the property that the 
Heisenberg algebra is the underlying algebra of the system. This con- 
struction can be generalized by considering formally the Hamiltonian 
as an operator acting in the g-deformed phase space. It implies the q- 
deformed Heisenberg algebra as the underlying algebra, xp — qpx = —i 
of the system and taking a class of transformations preserving this 
algebraic structure. It leads naturally to the Fock space formalism as 
an adequate language of our study. 

The harmonic oscillator as well as anharmonic, both classical and 
quantum plays a fundamental role in physics. A particular goal of 
present study is to find discrete systems those can be related to quantum- 
mechanical (an)harmonic oscillator via a quantum canonical transfor- 
mation. It will be considered two types of discrete systems: translation- 
covariant (uniform grid) (i) and dilatation-covariant (exponential grid) 
(ii). In what follows it will be used a Fock space formalism. In the 
next paper the case of Coulomb, Poschle- Teller and Morse potentials 
will be presented 0. 

Important general remark should be made in row. Generically, if 
two systems I and II are related by a gauge transformation and/or by 



^It also allows a possible modification of the Hamiltonian by inserting the pa- 
rameter q in appropriate places 
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a change of variables as equivalent, namely, 

l-ii{x) = g{y)^'^^l<n{y)g{y)\y=y(x) 
these systems will be treated as equivalent. 

1. Harmonic oscillator in Fock space 
The Hamiltonian of harmonic oscillator is defined by 

where u is the oscillator frequency. The configuration space of (1) is 
the whole real line, x G (—00, 00). The eigenf unctions and eigenvalues 
are given by 

^fc(x) = Hk{V^x)e-^'^ , Ek = uj{k + 1/2), A; = 0, 1, . . . (2) 

where Hk{y) is the fcth Hermite polynomial in standard notation. With- 
out a loss of generality we put all normalization constants in (2) equal 
to 1. The Hamiltonian (0) is Z2-invariant, x — x, which leads to two 
families of eigenstates: even and odd, or, in other words, symmetric 
and anti-symmetric with respect of reflection, correspondingly. This 
property is coded in parity of the Hermite polynomials and is revealed 
by the relation: 



H2n+piV^x) = uj2x'L'^ ~^\ujx^) , n = 0, 1, . . . (3) 

where L^\y) is the nth associated Laguerre polynomial in standard 
notation, and the parameter p = 0, 1 determines parity, P = (— 1)^. 
Hereafter we can call 

¥^\x) = xPe-'^'^ , (4) 

the ground state (the lowest energy state) of parity P. Thus, the 
formula (4) makes an unification of both possible values of parity and 
for the sake of simplicity we will call (4) the ground state eigenfunction 
without specifying parity. 

Let us make a gauge rotation of the Hamiltonian (1) taking as a gauge 
factor the ground state eigenfunction (4) and changing the variable x 
toy = ujx^, which incorporate the reflection symmetry. After dropping 
out non-essential constant term, (| + p), which cause a change the 
reference point for energy only, we get an operator 

h{y,dy) = - i¥^\x)r'n¥^\x) \y=^,2 = -2ydl + 2{y-p-\)dy 
^ ^ (5) 
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with the spectrum En = 2n, where = 0, 1, 2 . . . 0. Now the eiegen- 
value problem for the operator is defined on the half-hne y G [0, oo). 
The operator (5) simultaneously describes a family of eigenstates of 
positive parity if p = and a family of eigenstates of negative par- 
ity if p = 1. We will call (5) the algebraic form of the Hamiltonian 
of the harmonic oscillator. The word 'algebraic' refiects the fact that 
the operator (5) has a form of linear differential operator with polyno- 
mial coefficients and furthermore possesses infinitely-many polynomial 
eigenfunctions. The latter implies that any eigenfunction can be found 
by algebraic means by solving a system of linear algebraic equations. 

The algebraic form (5) admits a generalization of the original Hamil- 
tonian (1) we started with. If we assume that the parameter p can take 
any real value, p > —1/2, one can make an inverse gauge transforma- 
tion of the operator (5) back to the Hamiltonian form and we arrive 
at 

2yd'y + 2{y-p~^)dy 

,2 



= , (6) 



2"" ' 2" ' 2x2 

which is known in literature as Kratzer Hamiltonian. It is worth to 
mention that this Hamiltonian coincides with 2-body Calogero Hamil- 
tonian. Also this Hamiltonian appears as a Hamiltonian of the ra- 
dial motion of multidimensional spherical-symmetrical harmonic oscil- 
lator. Hereafter we will continue to call the system characterized by 
the Hamiltonian (6) the harmonic oscillator. 

The resulting Hamiltonian (6) is characterized by the eigenfunctions 



,2 



^lJk{x) = xPLi: -^'{ujx^)e~^^ , (7) 

which coincides with (2) at p = 0, 1. The spectrum (6) is still equidis- 
tant with energy gap uo and after appropriate shift of the reference point 
it coincides with the spectrum of the original harmonic oscillator (1). 
Thus, the deformation of (1) to (6) is isospectral which is, of course, 
well-known. 

In order to move ahead let us introduce a notion of the Fock space. 
It will be a natural formalism to study canonical transformations. Take 
two operators a and h obeying the commutation relation 

[a, h]= ah — ha = I, (8) 



^It is necessary to emphasize that the spectrum is defined by action of the op- 
erator 2ydy in (5) on monomial: {2ydy)y" = 2ny". It is the only operator of 
zero-grading in (5) 
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with the identity operator / on the r.h.s. - they span a three-dimensional 
Lie algebra which is called the Heisenberg algebra /is. By definition the 
universal enveloping algebra of /13 is the algebra of all normal-ordered 
polynomials in a, b: any monomial is taken to be of the form b^a"^ 
If, besides the polynomials, all entire functions in a, b are considered, 
then the extended universal enveloping algebra of the Heisenberg alge- 
bra appears or in other words, the extended Heisenberg- Weyl algebra. 
In the (extended) Heisenberg- Weyl algebra one can find the non-trivial 
embedding of the Heisenberg algebra: non-trivial elements obeying the 
commutation relations (^), whose can be treated as a certain type of 
quantum canonical transformations. We say that the (extended) Fock 
space, T is determined if we take the (extended) universal enveloping 
algebra of the Heisenberg algebra and attach to it the vacuum state 
|0 > defined as 

a|0 > = . (9) 

It is easy to check that the following statement holds. If the operators 
a, b obey the commutation relation (H), then the operators 



Jl = b'a - 


nb , 


Jl = ba- 


n 




2 ' 


J~ = a 





span the s/2-algebra with the commutation relations: 

[J°,J±] = ±J±, [J+,J-] = -2j\ 

where n G C For the realization (jT^) the quadratic Casimir operator 
is equal to 

C^2 = ^{J+J-}-« = -|(^ + l) , (11) 

where {, } denotes the anticommutator and is c— number. If n G Z+, 
then (|l^) possesses a finite-dimensional, irreducible representation in 
Fock space leaving invariant the linear space of polynomials in b acting 
on vacuum: 

Vnib) = (1,6,6^...,6")|0), (12) 

of dimension dimP^ = (n+1). It is evident that any operator which is a 
polynomial in generators preserves the space Vn{b) (and converse 

is also correct 0). Such an operator we call SI2- quasi- exactly- solvable 

■^Sometimes this is called the Heisenberg- Weyl algebra 
^For details and discussion see, for example, 
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operator. Thus, the most general s/2-quasi-exactly-solvable operator in 
the Fock space having a form of a polynomial in a of degree not higher 
than two is given by 

T2 = c^+J^J^ + c^qJI^J^ + Coo + Cq-J^J^ + c — J„ + 

c+J+ + CoJ° + c_4- + c, (13) 

where arbitrary c-numbers, or after the substitution ([ToD 

in explicit form 

T2(6,a) = -P^{b)a^ + P^{h)a + P2(&), (14) 

where the Pj{h) are polynomials of jth order with coefficients related 
to Cq,/3, Cq,, c and n. 

The spaces P„ possess a property that P„ C Vn+i for each G Z+ 
and form an infinite fiag (filtration) and 

Hence it is evident that any operator which is a polynomial in gen- 
erators J°'~ only preserves the fiag of P ||^. Such an operator we 
call SI2- exactly- solvable operator. Thus, the most general s/2-exactly- 
solvable operator in the Fock space having a form of a polynomial in a 
of degree not higher than two is given by 

E2 = Coo + Co-J°J' + c^-J',r + CoJ° + c_ J" + c , 

(15) 

where J^'° = Jg ' and c arbitrary c-numbers. After the 

substitution ([T0|) in explicit form 

^2(6, a) = -g2(&)a' + Qi{b)a + Qo , (16) 

where the Qj (b) are polynomials of jth order with coefficients related 

to Cq,^, Cq,, C. 

Now we can introduce a notion of the spectral problem in the Fock 
space. Let L{b, a) is an element of the Heisenberg-Weyl algebra. By 
definition, to solve the spectral problem for the operator L{b, a) is to 
find a set of the elements {</>(&)} in the Heisenberg-Weyl algebra and a 
corresponding set of parameters {A} for those the equation 

L{b,a)(f){b)\0 > = X(j){b)\0 > {*) 

is fulfilled. We will call {(f){b)} and {A} the eigenfunctions and eigen- 
values, respectively. Attempting to study the spectral problem (*) 
immediately leads to a delicate question about convergence of the op- 
erator series. In order to avoid possible difficulties we will restrict our 
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consideration by the cases when the polynomial in b eigenfunctions 
appear. 

2. Translation- covariant discretization 

Take as an example the s/2-exactly-solvable operator of the following 
form 

/i(i)(6,a) = -2J°J- + 2J° -2{p+-)J- = -2ba^ + 2{b - p - -)a , 

2 2^^^^ 

where p is a parameter and J^'° = J^'^ (see (10)). Simple analysis leads 
to a statement that the eigenfunctions of hf are the associated Laguerre 

polynomials of the argument b, Ln ^ (6) and their eigenvalues, En = 
2n, n = 0,l,2,... g. 

As the next step we consider two different realization of the Heisen- 
berg algebra in terms of differential and finite-difference operators. 
A traditional realization of (|^) appearing in text-books is the so-called 
coordinate- moment um represent at ion : 

d 

dy 

(see, for example, the book by Landau and Lifschitz p), where the 
operator b = y stands for the multiplication operator : bf{y) = yf{y)- 
In this case the vacuum is a constant and without a loss of generality 
we put |0 > = 1. However, there exists another realization of (^ in 
terms of finite-difference operators. It is a finite-difference analogue of 

my- 

a = V+, b = {y + a){l - 6V_) , (19) 

where 

fiy±S)-fiy) 

^±f{y) = ^ , 

(cf.[^]) is the finite-difference operator. It represents what can be called 
a 6— discretization of the derivative: V± dy, if 6 ^ and, thus, 
which can be called 5— derivative 0. In general, 6, a can be any com- 
plex numbers and T'±(— 5) = Vzp{6). It is necessary to emphasize 
that from physical point of view the operator b in ([T9|) is nothing but 



a = — = dy , b = y , {li 



^As in (5) the eigenvalues are defined by action of the operator (26a) in (13) on 
monomials: (26a)&" = 2n6" + 2&"+^a. The second term disappears after action on 
the vacuum 

^This finite-difference operator is also known in mathematics literature as the 
Norlund derivative (see, for example, Q]) while we prefer to use a name S -derivative 
in order to distinguish from q-derivative (see below) - another type of discretization 
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the canonical-conjugate to a discrete momentum operator defined by 
5— derivative. So, (0) represents two-parametric family of quantum 
canonical transformations. For a = and 6 —>■ 0, the formulas (|T9|) 
become (|18|). It allows to interpret (^) as the continuous deformation 
of (0)- 

A remarkable property of this realization is that the vacuum can 
be taken as a constant [] and thus, without loss of generality, it can 
be placed as |0 > = 1 for both cases (|T8D-(|l9|). Realization ( ]T9| ) is 
translation-covariant: under a linear shift of variable y ^ y + B the 
functional form of (|T^ is preserved. 

Substitution of (|I^) into (|T^ leads to the operator (5) - the algebraic 
form of the Hamiltonian of the harmonic oscillator. Thus, the operator 
([17|) can be called the the Hamiltonian of the harmonic oscillator in the 
Fock space. It is evident that the procedure of realization of the Heisen- 
berg generators a, b by concrete operators (differential, finite-difference, 
discrete) provided that the vacuum remains unchanged leaves any poly- 
nomial operator in a, b isospectral. 

Now let us study another 'face' of harmonic oscillator by substituting 
the realization of the Heisenberg algebra by finite-difference operators 
into (0). It results to 

hf\y, Di) = -hy + a + 6{p + + 2(1 + ^){y + a)D_ . 

6 2 6 ^20) 

In this realization the corresponding spectral problem (*) can be writ- 
ten as 

-^[y + » + ^ip + 1)] (t^iy + ^) + ^ [(1 + + a + P + ^] </>(2/) 



- ^(1 + \){y + - 5) = Ecl>{y) . (21) 

It is worth to note that the Askey condition that the sum of the co- 
efficients in front of unknown functions should be zero is fulfilled at 
E = It is equivalent to a statement that the equation (|2T| ) possess 
an eigenfunction which is a constant. 

In general, the operator hs{y,D±) is a non-local, three-point, finite- 
difference and translation-covariant operator. This operator is canonically- 
equivalent to the harmonic oscillator. Pictorially it is illustrated by 
Fig.l. 



^Any 5— periodic function can be chosen as a vacuum for dl^) as well 
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(t){y-5) 4>{y) 4>{y + 6) 

Fig. 1. Graphical representation of the operator (|20| ) 

So, in the y— space we have the uniform grid, hnear lattice. However, 
from the point of view of the original configuration space, x— space 
where the harmonic oscillator (1) is defined we have a square lattice. 

In order to find the eigenfunctions of (|20| ) a certain trick can be used. 
One can easily show that the following equality holds 

[{y + a)e-''yri = {y + af^h , 

where {y + = [y + a){y + a — 5) . . . {y + a — nS) is a so-called 

quasi-monomial or generalized monomial and / is the identity opera- 
tor. Then making use this relation it is not difficult to prove that the 
eigenfunctions of (^Of ) remain polynomials in y and, furthermore, they 
can be found explicitly in rather elegant way 

n 

My) = Lt''\y,S) = ^a?-^\y + , (22) 

e=o 

where the coefficients in the expansion of the Laguerre poly- 

nomials, Ln ^\y) = X]"=o'^^^ ^^y^- '^^^^ these polynomials the 
modified associated Laguerre polynomials. Simultaneously, the eigen- 
values in the equation ([21| ) remain equal to (2n),n = 0,1,2... and 
they are the same as the eigenvalues of the harmonic oscillator prob- 
lem (1) as well as (6) and (13). Thus, one can say that the operator 
( |20| ) defines a finite- difference or 6— discrete algebraic form of harmonic 
oscillator Hamiltonian. Without loss of generality one can place a = 
in (^, (©, (g2D. 

There exists two non-trivial particular cases of (|20D. The first case 
corresponds to the spacing 6 = —1. It leads to a disappearance of the 
term proportional to D_ and, thus, the operator ( PP] ) becomes two- 
point operator (!) 

h^s\y,D±) = 2[y + a- {p + ^)]D+\s=., = 

2[y + a-{p+^)]D_\s=i . (23) 

It is worth noting that the spectrum in this case is defined by the 
term {2yD^\s=-i) stemming from the negative grading term (26a^) in 
(p!7|). Breaking a condition of performing the canonical transformations 
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we re-insert the parameter 6 in (^3[) in naive, straightforward manner: 
D-^-\s=-l D^{6). It resuhs to the operator 

hf\y,D±) = 2[y + a-{p+^)]D+ , (24) 



which is isospectral to (5), (|T^ for any 6. In limit 5 — it leads to 
the first-order differential operator 

hf\y,dy) = 2[y + a-ip+^)]dy. (25) 



Although, the operators (PHj), (^) are isospectral to various forms of 
the harmonic oscillator (1), (6), those operators are not related to each 
other by a quantum canonical transformation and therefore are not 
canonically equivalent. Above-mentioned isospectral transition from 
the second-order differential operator (5) to the first order differential 
operator (|25|) or from finite-difference one (^OD to (^4]) reminds cele- 
brated Bargmann transformation (see for example p|). It is necessary 



to emphasized that the operator (|2^) has infinitely-many polynomial 



eigenfunctions, however, unlike the operators (13), (16), the operator 
(20) is not s/2-exactly-solvable (!). Also it can not be rewritten in 
terms of the generators a, b of the Heisenberg algebra (8) and does not 
belong to the Fock space (see discussion below). 

Another important particular case occurs if the spacing 6 = —2. In 
this case the operator (pO]) again becomes two-point one 



(26) 



-^[y + a-2ip+ i)]0(y -2) + ^{y + a)<P{y + 2) = E(f){y) 

where a new spectral parameter E = E + a+ p+^. 

Remark. The function (f){y) in the r.h.s. of ( pij ) can be 
replaced by (j){y + 6) or (f){y — 6), or by a linear combination of 
0(?/±5), 4>{y). It does not change the statement that (plf ) has 
infinitely-many polynomial eigenfunctions. This procedure 
preserves isospectral property. However, such changes of the 
r.h.s. lead to a replacement of the original standard spectral 
problem (*) by a 'generalized' spectral problem. In this case 
the r.h.s. contains an operator other than the operator of 
multiplication on a function. In general, a physical relevance 
of such right-hand-sides is unclear. It is worth to note that it 
reminds the Sturm representation of the Coulomb problem 
when the energy is kept fixed but a set of discrete electric 
charges for which a corresponding excited state energy is 
equal to this energy is studied. In this formulation the r.h.s. 
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of the Coulomb problem contains the Coulombic interaction 
potential as the weight factor. 

A natural question can be posed about the most general second-order 
linear differential operator, which (i) is isospectral to the harmonic os- 
cillator (1),(6); (ii) has infinitely-many polynomial eigenf unctions and 
(iii) related to (5) by a canonical transformation. Following the Theo- 
rem one can derive the most general operator with above properties 

hs,g{y,dy) = -2{AJ^ + BJ-)J- + 2J^ -2CJ- = 

-2iAy + B)d'y + 2iy-C)dy (27) 
where A,B,C are arbitrary constants and the generators (|TDp are re- 



alized by differential operators (JTBI). However, by a linear change 



of variable, y — >• ay + P with appropriate a ^ 0, P the operator 
(|27|) can be transformed into (5). It reflects a fact that the linear 
space of polynomials 'P„(x) is invariant under a linear transformation: 
X 7iX + 72,7i 7^ 0. Thus, without loss of generality we can put 
A = a = 1 and also C = p + 1/2. The eigenfunctions of (^) remain 

the Laguerre polynomials but of a shifted argument, Ln {y + (3). It 
leads to a statement that among the second-order differential opera- 
tors there exists no non-trivial isospectral deformation of the harmonic 
oscillator potential preserving polynomiality of the eigenfunctions. 

The operator ( ^Tf ) can be rewritten in the Fock space formalism by 
using (|18D 



hg{h,a) = -2{b + B)a^ + 2{b - C)a . (28) 

It is evident that the operator (^) is the most general second order 
polynomial in a, which is isospectral to (O) and also preserves the 
space of polynomials (|12|). By substitution (|19|) into the operator (pSf) 
it becomes a finite-difference operator 

hs,g{y,D±) = -2BDl -hy + a + SC]D+ + 2(1 + hiy + a)D^ 

^ ^ (29) 



(cf. (pUD). Thus, in the realization by finite-difference operators the 
corresponding spectral problem (*) has the form 

Z^^{y + 25)-^[y + a-2B + 5C]<P{y + 5) 
+ ^[{l + -^){y + a)-j + C]<j>{y)--^{l + -^){y + a)<j>{y-6) = E<p{y) 



(30) 
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and is characterized by existence of infinitely family of eigenfunctions 
given by polynomials in y. 

The operator hs{yi D±) now becomes the four-point finite-difference 
operator, see Fig. 2 . 



0(1/ -5) 0(y) <P{y + 5) 0(1/ + 25) 

Fig. 2. Graphical representation of the operator (pQ]). 

It is quite surprising that a simple transformation like linear shift 
of variable y in differential operator ( pO|) (which leads to nothing non- 
trivial, see above) leads to occurrence of the extra point in the isospec- 
tral finite-difference counterpart changing a type of its non-locality. 
Nevertheless, these operators remain to be canonically equivalent being 
related each other by a canonical transformation. So, by a canonical 
transformation one can change a non-local nature of finite-difference 
operators. 

It is quite interesting to abandon the condition of canonical equiv- 
alence and pose a question about the most general differential (finite- 
difference) operators isospectral to the harmonic oscillator and possess- 
ing the infinitely-many polynomial eigenfunctions. One can easily show 
that ignoring the condition of canonical equivalence leads to nothing 
new for differential operators. However, for the case of finite-difference 
operators much wider class of operators occurs than ( pOD or (pOj). As 
a natural constraint we have to impose a condition of maximal num- 
ber of points in the finite-difference operators we search for. Let us 
consider three-point operators (see Fig.l). A simple analysis leads to 
a statement that the three-parametric operator 

hs{y, D^) = 2[(1 - A)y + C+)]D+ + 2[Ay + C^)D_ , (31) 

where A, C± are parameters, is the most general three-point finite- 
difference operator with infinitely-many polynomial eigenfunctions, which 
is isospectral to the harmonic oscillator. Of course, the operators (16) 
and (18) are particular cases of (|3lD . The form ( pi]) can be understood 
as a consequence of the Theorem which states that any finite- difference 
operator hs{y,D±) preserving the infinite flag of polynomials 

VM = (l,l/,l/^...,l/") , (32) 

should have a representation in terms of the operators: 



J± = , 4 = yD± 



(33) 
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As a remark we should note that the generators (pB]) do not span an 
algebra closed with respect to commutators. However, if we consider a 
linear space spanned by (|33|), then the Heisenberg algebra ( [T8| ) appears 
as its subspace. 

It is worth to mention that the canonical discretization of the har- 
monic oscillator can be made directly in the configuration a;— space 
where the original harmonic oscillator (1) is defined. If we make a 
gauge transformation (5) with (4) at p = as the gauge factor but 
without the change of variable, then 

is another algebraic form of the harmonic oscillator (cf.(5)). The opera- 
tor (1^) has the Hermite polynomials, Hk{\/ujx), as the eigenf unctions. 
Now we can make 5— discretization firstly rewriting (0) in the Fock 
space formalism 

h^^\h,a) = -^a^ + ba = -^J'J- + J\ (35) 

(cf.(|17[)) and then realizing a,b by finite-difference operators ( plQ]) 

hf\x,D±) = -^Dl + {x + a)D_ . (36) 
2iu 

This operator is canonically-equivalent to the harmonic oscillator, it is 
defined on the uniform linear lattice in x— space and has the eigenvalues 
Ek = k, k = 0,1,2 ... . However, it does not possess the symmetry 
X —X unlike the 6— discretized, canonically-equivalent operator ( pO]) 



in y— space. The operator (|36| ) is non-local, four-point operator with 
eigenf unctions 

k 

= (a; + «)^'^ (37) 

£=0 

where z^^^ is quasi-monomial and ae are the coefficients in the expansion 
of the Hermite polynomials, i7fc(x) = Yle=o^^^^- These eigenfunctions 
are closely related to the polynomials which can be called the modified 
Hermite polynomials 

k 

^fc(x,5) = 5^a,xW , (38) 
e=o 

and, thus, 

0fc(x) = Hk{^/uJ{x + a),^/uJ6) . 
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Summarizing, in this Section we presented two 6— discretized oper- 
ators which are canonically-equivalent to the harmonic oscillator, (pO]) 
and (|36|), which are defined on the uniform grid in y— and x— spaces, 
respectively. In spite of the fact that they are represented by differ- 
ent elements of the Heisenberg-Weyl algebra, they have infinitely-many 
polynomial eigenfunctions and are isospectral to the original harmonic 
oscillator (1). 

3. Dilatation- covariant discretization 

Apart from the translation-covariant discretization given by 5-derivative 
(see (0)) there exists a dilatation-covariant discretization based on q- 
derivative, T>q, which is also called the Jackson symbol 

where g is a complex number. This Section will be devoted to a brief 
discussion of the dilatation-covariant discretization or g— discretization. 
First of all by following the above-mentioned philosophy a natural ques- 
tion can be posed about an existence of a quantum canonical conjugate 
to g-derivative - the operator which obeys together with Vg the com- 
mutation relations (H). Up to our knowledge a definite answer is not 
found so far and very likely such an object does not exist in terms 
of well-defined operators. However, it is well known that the deriv- 
ative T>g appears naturally in connection to a realization of quantum 
algebras in action on functions in one and several variables. Thus, it 
looks reasonable to explore a generalization of the underlying Heisen- 
berg algebra described above to the case of the g— deformed (quantum) 
Heisenberg algebra. Hence, we will study 'deformed' quantum systems 
possessing a g— deformed (quantum) Heisenberg algebra as a hidden 
algebra instead of the standard Heisenberg algebra (H) (see discussion 
in Introduction). 

In order to proceed let us ask first what would happen if in the 
expressions ([I^), ( P^ ) the operators a,b are not the generators of the 
Heisenberg algebra (|]) but the generators of the g— deformed Heisen- 
berg algebra 

[a, b]g = ah — qha = 1 , (39) 

where g is a parameter. Following the Theorem proved in ||^, one 
can demonstrate that within the g— deformed Fock space built on the 
g— deformed Heisenberg algebra (|39D there exists the fiag V of linear 
spaces of polynomials in b (see (12)), which is preserved by the opera- 
tors (0), (|28|). By a simple calculation one can find the eigenvalues of 
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the operators (p]?]), ( p8D in the spectral problem (*) 

Ei") = 2M , n = 0,l,... , (40) 

where 



is a so-called g— number and {n} — > ra, if g — > 1. It is evident that if 
the parameter g is a real number the spectra of (|T^, (^Sj) are real. 

The algebra (^) has a realization in terms of g— derivative and the 
operator of multiplication (see, for example, @]) 

a = V„ b = y , (41) 

with the same g as in (^). This realization has a property that the 
vacuum remains the same as in the cases ([T8| ) - ([T9|) and without loss of 
generality it can be set as |0 > = 1. Now we can substitute ( PD in 
(17) and the following operator emerges 

h^\y,V^) = -2j0j- + 2j0-2(p+l)J- = 

-2yV'^ + 2{y-p~^)V,, (42) 

where the generators J° = ba, J~ = a have the same functional form 
as in (p!OD but obey the g— deformed commutation relation 

g 

forming the g— deformed Borel subalgebra b{2)q of the g— deformed al- 
gebra sl{2)q (for discussion see 0). In this case the operators (|17D, 
(|28|) are the sl{2) q-exactly-solvahle operators. Moreover, the operator 
(p!7|) (as well as (pS])) can be called the q— deformed harmonic oscillator 
Hamiltonian in Fock space possessing sl{2)q hidden algebra. 

The operator hq{y, Dq) is a non-local, three-point, discrete, dilatation- 
covariant operator defined on exponential lattice. It is illustrated by 
Fig.3. 



(t>iy) 0(^1/) 0(?^i/) 



Fig. 3. Graphical representation of the operator (|4^ ) 



The operator (^21) can be called the algebraic form of the Hamiltonian 
of the q~discretized harmonic oscillator. 



15 



The spectral problem for the operator ( ^21) has a form 



yq{q 



2 + q + - 2pq{l - q) 1 
q{q-lf y 



<P{<iy)- 



l + q- 2p(l - q) 

y{q-i? 

or, the r.h.s. can be taken as 
or as 



+ 



l-q 



E^^'^<P{qy) . 



E^'^Uiq'y) 



(43) 



(44) 



(45) 



Usually, the spectral problem for g— discrete operators is defined with 
(|4^) as the r.h.s. (see, for example, [|10|)- K assumes that the middle 
point in Fig. 3 remains fixed under dilatation. Introducing the new 
variable y = qy, it can be seen explicitly. 

If in the case (^3]) the eigenvalues are given by ( ^OD while for (^4|), 
(H5|) the eigenvalues are equal to 



i^W = _2{-n} 
E^^^ = 2q-'^n} 



n 



n 



0,1, 



0,1, 



(46) 
(47) 



correspondingly. In the limit g — > 1 all three expressions coincide 
corresponding to the original harmonic oscillator spectrum. The spec- 
tral problems (^)-(|45|) can be considered as a possible definition of 
a g— deformed harmonic oscillator. In the literature it is well-known 
many other definitions of the g— deformed harmonic oscillator (see for 
example. 



11 , 12 



and references therein, |]13l) 0. Such a situation re- 
flects an ambiguity of making a g— deformation | as well as absence 
of clear physical criteria, which can remove or reduce this ambigu- 
ity. For instance, in the literature it is exploited three different types 
of the g— Laguerre polynomials (see, for example, an excellent review 
||14||), but it is not clear why other possible g— deformations of Laguerre 



polynomials are not studied. 



^Usually, these deformations are done by a direct discretization of the origi- 
nal Hamiltonian (1). Most of all are based on a discretization of the Infeld-Hall 
factorization representation of (1) 

^For example, any term in non-deformed expression can be modified by mul- 
tipliers of the type and even some extra terms can be added with vanishing 
coefficients in the limit q — > 1 like (1 — g)**, 6 > 



16 



Isospectrality of (17) and (28) is preserved by the g— deformation. 
Substitution of (|4l|) in (p8|) gives a slight modification of the expressions 
(|i2|). Unhke translation-covariant case it does not lead to a change 
of the nature of non-locality changing the number of points in the 
operator (^2|) as it is happened for the operators (^) and ([29|). The 
g— deformation of another algebraic form of the harmonic oscillator 
(|3^ ) in Fock space 

hf\h,a) = -1-J-J- + J\ (48) 



hf\x,V,) = -—pi + xV,, (49) 



(cf.(^5])) takes in terms of the g— derivative the following form 

1 

2uj' 

(cf.(^2])). It should be mentioned that the operators (^2]) and (|^) are 
defined on the essentially different lattices, which are exponential in 
X— and ?/— variables, respectively. 

Similar to what was done for canonical transformations it seems nat- 
ural to introduce a notion of g— deformed canonical transformations, 
when two g— deformed systems are g— canonically equivalent if they can 
be connected through the g— deformed canonical transformation (see 
below). However, we were unable to find well-defined, non-trivial real- 
ization of the algebra ( ^9]) other than (^) which, for instance, would 
be similar to the realization f|20|) for flSl) (see discussion in 01) . 



4. Anharmonic oscillator (perturbation theory) 

Anharmonic oscillator is one of the most important non-exactly- 
solvable problems of quantum mechanics. One of the simplest and 
the most popular examples is given by the Hamiltonian 

= + Y^' + , n = 2, 3, 4, . . . (50) 

where g is a. coupling constantf^ Making first the gauge rotation of 
(|50|) with the gauge factor @) we get 

UJ 

= -2ydl + 2{y-p- i)9, + Ay" , (51) 

(cf.(^), where p = 0, 1 and has a meaning of parity and A = guj~^'^^^\ 
By the reasons which will be clear later the operator (pTI) can be called 



"'^''So far, the most comprehensive study of the anharmonic oscillator ( pO| ) was 
carried out by Bender- Wu in the classical paper [|l5| (see also [|6|) 
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algebraic form of the anharmonic oscillator Hamiltonian. Replacing in 
(|51|) the derivative and the coordinate by the elements of the Heisen- 
berg algebra (|^), — a, ?/ — 6 we will arrive at the element of the 
Heisenberg-Weyl algebra 

h^''^°\h,a) = -2ba^ + 2{b - p - + A6" , (52) 

which can be called the anharmonic oscillator Hamiltonian in the Fock 
space. Now we can fix n = 2 and we will consider this particular case 
of the quartic anharmonic oscillator as the major example in further 
consideration. 

Let us find a finite-difference operator which is canonically equivalent 
to dSPD- In order to do it we simply substitute the realization (|TP|) of 
the Heisenberg algebra to the operator (^) 

ht'"\y,D^) = -'^^[y + a + S{p + ^)]D^ + 2[{1 + -^){y + a) - 

^(y + af^]D^ + ^(y + + ^iy + af\ (53) 

(see Fig. 4). It is quite amazing that the perturbation leads solely to 
an addition of one more point (marked by o) to the harmonic oscillator 
operator. A model characterized by the operator (^) can be called a 
finite- difference anharmonic oscillator. 



<P{y-2S) ^(y-S) 0(y) <P{y + S) 

Fig. 4. Graphical representation of the problem (^) 

It is non-local, four-point, finite-difference operator (for comparison see 
Fig.l and Fig. 2 with S —S). 

The corresponding spectral problem (*) looks as follows 

-^[y + (^ + ^ip + lWy + ^)+ l[i'^ + l)y + (^ + p + l]4>{y) - 

hl + h(y + a)<P{y-5) + A(y + a)(2)0(y-25) = E<P{y) . 

^ ^ (54) 

Their eigenvalues coincides with those of the anharmonic oscillator 
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We have to emphasize that a presence of the anharmonic term changes 
the nature of non-locahty of the harmonic oscillator leading to appear- 
ance an extra point. In particular case 6 = —1 the number of points is 
reduced to three, however, the lattice becomes non-uniform. 

In similar way one can construct a g— deformed anharmonic oscillator 
taking the operator (^) as an element of the g— Fock space. Substi- 
tuting the realization ( PT| ) of the g— deformed Heisenberg algebra into 
we get 

= -2yVl + 2{y-p~^)V, + Xy\ (55) 

(cf. (0)) which three-point, non-local operator (see, for example, 
Fig.3). 

The spectral problem (*) for the operator (|55D has a form 



yq{q - l)^ yq{q - 1) 

2-(2p+l)(l-,) + 2,(l-.)-A,^fa-l)^ = E>'>Hy) . 

yig - ir (56) 

(cf.(^2])). Certainly, the r.h.s. in (|56D can vary being equal to either 
(|4^), or (^5]). It reflects a fact that in the case of the g— Fock space the 
spectral problem (*) can be modified in one way 

L(6,a)0(6)|O > = A0(g6)|O > , (**) 

or another 



L(6,a)0(6)|O > = A0(g'6)|O > 



* * * 



If we introduce the multiplication operator Tgf{x) = f{qx) the prob- 
lems (**) or (***) correspond to an appearance of non-trivial operator 
weight factors in r.h.s. 

L{b,a)(j){b)\0 > = XTq(f){b)\0 > , 

or 

L(6,a)0(6)|O>= Ar20(6)|O>, 

respectively. 



It can be easily shown that the operator (|52D describing the an- 
harmonic oscillator does not belong to the class of the s/ (2)— exactly- 
solvable operators. Hence their eigenfunctions are not polynomials. 



However, it was proved in [|T^ that in framework of some perturbative 
approach (see below) as a consequence of the fact that the perturba- 
tion 6^ G V2{b) (see (|1^) the perturbation theory in powers of the 
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parameter A is algebraic one: any correction to any eigenfunction is a 
finite-order polynomial and hence can be found by algebraic means. In 
fact, such a perturbative approach provides a certain regular way to 
define an object in the Heisenberg-Weyl algebra which can be called the 
Hamiltonian of the quartic anharmonic oscillator. From practical point 
of view such a perturbation theory has very important feature: once 
being developed for the operator ( |52D it gives a unique possibility to 
construct simultaneously (!) the perturbation theory for the operators 



Let us approach to a construction of the above-mentioned pertur- 
bation theory. Following a standard prescription we take the spectral 
problem (*) with 

L{b, a) = ho + Xhi 

and develop a perturbation theory in powers of A searching for correc- 
tions in a form 

</. = 5];A>„ , E = ^A"E„, . (57) 

Collecting the terms of the order A" it is easy to derive an equation for 
the nth correction 

n 

{ho - -Eo)0n = X] Ei(pn-i - hi(j)n-l ■ (58) 
i=l 

A remarkable feature of this form of perturbation theory a possibility 
to study a single state separately, without touching other states as it 
was the case of the Rayleigh-Schroedinger form of perturbation theory. 

Now we take as the unperturbed operator Hq the Hamiltonian of the 
harmonic oscillator in the Fock space ([T7|) and consider as the perturba- 
tion hi = b"^. As it was mentioned already, in this case the nth correc- 
tion (pn to eigenfunction should be polynomials in b. As an instructive 
example let us study the ground state of anharmonic oscillator. The 
ground state of the unperturbed problem (O) is characterized by 



00 = 1 , ^0 = 1 • (59) 

By solving the equation (|58D after simple calculations we can get the 
explicit form for the first several corrections, for example, 

b' , (3 + 2p) ^ (l + 2p)(3 + 2p) 
-'^^ = 2{2} + ' = 4 ' 



and 

h = \ 



^ 2g^ + 5g + 6 + 2p(g + 2) ^3 ^ 



{2}{4} 2{2}{3} 
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+ 12g + 15 + 8p{q + 2) + V ^ (3 + 2p)[g^ + 3g + 3 + 2p(g + 1)] 
4{2} 2 

= _ (l + 2p)(3 + 2,) _ ^g^^ 

o 

Without any difficulties one can find several next corrections. However, 
it becomes evident very quickly that the complexity of calculations is 
growing very fast with a number of correction. Using different real- 
izations of the (g) — Heisenberg algebra one can calculate perturbative 
corrections to the various forms of anharmonically-perturbed harmonic 
oscillator (differential, finite-difference, discrete). 



I. q 



1, b 



This case corresponds to the coordinate-momentum realization of 
the Heisenberg algebra (|r^) and vacuum definition |0 >= 1. It 
leads to a standard anharmonic oscillator (BDI) and the corrections 



are: 



and 



(3 + 2p) 

4 4 ^ 



y , ii + 6p^3 



+ 



12 



31 + 24p + 4p2 



y' + 



(3 + 2p)(7 + 4p) 



:i + 2p)(3 + 2p)(7 + 4p) 
8 



(61) 



where Ei is the same as in (|60|). 

Given form of the perturbation theory coincides to the so-called 
'F— functions method' developed by Dalgarno (see discussion in 
p6|), which in fact was realized for the case of the ahnarmonic 
oscillator ( pOD in |T^ . It is easy to check that the corrections ( |6lD 
coincide to those calculated in text-books (see for example [0). 



II. g = 1, 6= (y + a){l - SV_] 



see 



This case corresponds to the translation-covariant discretization 
and perturbed finite-difference harmonic oscillator (^) (see (|53|) ): 

y(^) , (3 + 2p) , 
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and 



8 12 ^ 



31 + 24p + V^^^2) 



(3 + 2p)(7 + 4p) , 



(62) 



where i?i is the same as in (|60D and E2 is the same as in (|6TD. 
Here = (2/ + = + «)(?/ + a - 5) ...(?/ + a - is 

quasi-monomiaL 

III. q^l,b = y. 

This case corresponds to the perturbed g— harmonic oscillator 
(|55|)-(|56|) and corrections are equal to 

■2 (3 + 2p) ^ (l + 2p)(3 + 2p) 



y 



+ 



-y , El 



4 



and 



{2}{4} 



2/^ + 



2g^ + 5g + 6 + 2p(g + 2) 
2{2}{3} 



r + 



4g2 + 12g + 15 + 8p{q + 2) + 4^^ (3 + 2p)[g2 + 3g + 3 + 2p{q + 1)] 



4{2} 



;i + 2p)(3 + 2p) 



[g2 + 3g + 3 + 2p(g + 1)] . (63) 



It is quite interesting to see how above-mentioned results will be 
modified if instead of the spectral problem (*) the spectral problem 
(**) is considered. The equation (p3) for the nth correction becomes 



{ho - EoTg)(f)n = ^ EiTg(j)n-^i - hi(j)n-l 



(64) 



i=l 



(cf. (|58D). The ground state of the unperturbed problem (|59D is un- 
changed. The first correction is also unchanged while the second cor- 
rection now takes a modified form 



{2}{4} 



^4 ^ 2g^ + 5g + 6 + 2p(g + 2) ^3 ^ 



2{2}{3} 



4g2 + I2q + 15 + 8p{q + 2) + Ap\^ (3 + 2p)[q^ + 3g + 3 + 2p{q + 1)] . 



4{2} 
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^ _ (l + 2p)(3 + 2p) ^ 8p(7-3g^) + 12^^(1 -g^)] . 

64 (65) 

5. Anharmonic oscillator (quasi-exactly-solvable model) 

Among one-dimensional Schroedinger equations there is some class 
of problems possessing a certain outstanding property - first several 
eigenstates can be found explicitly, by algebraic means. Such problems 
are called quasi-exactly-solvable [0. Among ten known families of 
one-dimensional quasi-exactly-solvable potentials |T3 there exists one, 



which can be treated as an anharmonic oscillator and its Hamiltonian 
is 

n^'-'^ = -^1^ + [y - (2^ + ^ + P^aW + gu^x' + ^^x' ■ 

^ (66) 

Here is a coupling constant and x G (— C)0, oo). Parameter p = 0, 1 
has meaning of parity. The first in + 1) eigenfunctions of parity p (but 
not others) are of the form 

^(f)(x) = x^'P„(cux2)e-""'/2-5xV4 

where Pn{y) is a polynomial of degree n. 

Making first the gauge rotation of (^Up with the gauge factor ( pTj ) at 
ri = we get 



to 



X) \y=uJX^ 



= -2ydl + 2{\y^ + y-p-l)dy - 2\ny , (68) 

(cf. (^,(^11)), where A = goj^'^ and constant terms are dropped out. 
The spectral problem for ( |68D is defined on the half-line, y G [0, cxo). 
The first (n + 1) eigenfunctions of (|68|) are some polynomials of the 
degree n, Pniy) (cf. (|67|) ) possessing /c = 0, 1, . . . , n real zeroes inside of 
the interval y G [0, oo). 

Replacing in ( |68D the derivative and the coordinate by the elements 
of the Heisenberg algebra (^: dy ^ a,y b we arrive at the element 
of the Heisenberg- Weyl algebra 

/i(9^^)(6,a) = -2ba^ + 2{\b^ + b-p-^)a - 2\nb , (69) 

(cf . (0) , ( p8D , (p2D ) . As in previous consideration the operator h^''^^^ {b,a) 
can be treated as an element of the Fock space as well as an element of 
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the g— Fock space. This operator can be called the Fock space Hamil- 
tonian of the anharmonic quasi- exactly-solvable oscillator. The opera- 
tor h^i^'^b, a) in the Fock space is s/2— quasi-exactly-solvable operator. 
It can be rewritten in terms of the generators of s/2— algebra (p^) 

h(^^^\b, a) = -2 J° J- + 2AJ+ + 2 jO - {n + 1 + 2p)J- , 

(70) 

where as always the constant terms are dropped off. 

However, it is easy to check that in the g— Fock space the operator 
a) is not 5/2^— quasi-exactly-solvable operator, since it can not 
be rewritten in terms of the generators of 5/2^— algebra Q. This operator 
needs a slight modification of the last term in order to become the 
s/2g— quasi-exactly-solvable: 

h^'J'^'\b,a) = -2ba^ + 2{\b^ + b - p - -)a - 2\{n}b = 

2 (71) 

-2 J° + 2A J+ + 2 J° - f 2Mi^ + l + 2p] , 



{2n + 2} 

(cf. ([69|)), where {n} is the g— number (see (31)) and the constant 
terms in the second expression are dropped off. 

The first {n + 1) eigenfunctions of (71) are some polynomials in b of 



degree n, Pnip). If in the case n = the first eigenfunctions for the 
spectral problems (*) and (**) coincide and equal to: 

= 1 , = , 



while for the case of n = 1 they are already different. Namely, for the 
spectral problem (*), 



4^) 



It Vl + 4A(l + 2p) 



4A 



^^In the case of the sZ2(j— algebra the major modification of ( |l0| ) comes for the 
positive-root generator 

= b^a — nb = b^a — {n}b , 

while the Cartan generator 

•^'^^^ 2 " {2n + 2} ' 

where {n} is the g— number (see (31)) and remains unchanged; for details, see 

13 
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where sub-indices (+) and (— ) are assigned to the ground and the 
first excited state, respectively. As for the spectral problem (**) the 
formulas (1721) are modified 



,(1) , It Vl + 4Ag(l + 2p) 
=b 4A ' 

^(1) ^ 1 T yi + 4Ag(l + 2p) _ ^^^^ 

Of course, when q = 1 the formulas ([72|) and (^) coincide. 

Now we can find a finite-difference operator canonically equivalent 
to ([681). In order to do it we put g = 1 in ( p8| ) and substitute the 
realization (^) of the Heisenberg algebra into the operator (p9| ) 

~[y + a + Sip+l)]D+ +2[{l + h{y + a) + Xiy + a)iy + a + S{n~l))]D_ 
old 

-2X6{y + a)^^^D_D_ - 2Xn{y + a) , (74) 

(cf. (p3|) ). It turns out that this operator is four-point finite-difference 
operator (see Fig. 4). It is quite amazing that independently on n the 
perturbation 2Xb{ba — n) to the harmonic oscillator operator (see Fig.l) 
leads solely to an appearance of one extra point (marked by o in Fig. 4) 
similarly to (|53|). A model characterized by the operator (^) can be 
called a finite- difference quasi- exactly-solvable anharmonic oscillator. 

Without loss of generality we place a = in ( [74D and the corre- 
sponding spectral problem (*) for (173) looks like 



-\[\ + p + \]<P{y + ^)+ -^[{^ + \)y + p + \\ <P{y) 



[l + ^) - X{y - 6{n + l))]<P{y - 6) 

- 2X^<P{y-26) = E4>{y) . (75) 



It is natural to assume that the spectral problem (|75D is defined in 
y G [0, oo). Their first (n+1) eigenvalues coincide with those of the an- 
harmonic oscillator Hamiltonian (^) as well as the operators (|6^) , (|69D . 
Their eigenfunctions remain polynomials but modified - each monomial 
should be replaced by quasi- monomial, y'' y^^\ For instance, the 
formulas ([721)- (ff^) become 

,(1) _ It Vl + 4A(l + 2p) 

^± - y ' 
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4" = lTv/l + 4A(l + 2p) 

It is worth to note that for physically relevant parameters g > and 
q > 0the eigenfunctions ([72D, (|76D in domain y >0 have no nodes for 
ground state and have the only one node for the first excited state like 
an analogue of the oscillation (Sturm) theorem holds (see, for example, 
PI). Likely, an agreement with Sturm theorem will hold for higher 
excited states but we were unable to prove it in full generality. 

6. Conclusion 

We introduced a simple-minded notion of canonical equivalence in 
quantum mechanics. Then we showed that for canonically-equivalent 
systems the spectra remain the same and constructed discrete systems, 
which are canonically-equivalent to the harmonic and a certain type of 
anharmonic oscillators - the quartic oscillator and the quasi-exactly- 
solvable oscillator. In general, these discrete systems are defined on 
non-linear lattices. 

We restricted our studies to the operators in the algebraic form pos- 
sessing the polynomial eigenfunctions. An important question to ask is 
what would happen if the operator in hand possesses an algebraic form 
but non-polynomial eigenfunctions. Perhaps, the most explicit instruc- 
tive example is given by the original Hamiltonian (1) of the harmonic 
oscillator. Following our philosophy it can be rewritten as the element 
of the Fock space 

n = -la^ + ^b\ (77) 
where its ground state eigenfunction can be written formally as 



oo 



Mb)\0>= E^^'"|0> • (78) 

n=0 

Substitution of (19) with |0 >= 1 into (78) leads to an infinite series 
of the form 



n=0 



which has zero radius of convergence in x for any 5 7^ 1^. Although 
the Taylor expansion in powers of x for the original eigenfunction (4) 



^^In this case q has the meaning the parameter of dilatation 
^•^This type of expansion on the basic set of quasi-monomials is known in htera- 
ture as the Newton series. For discussion see, for example, 
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at p = had infinite radius of convergence. So, the radius of con- 
vergence has delta-function behaviour. Perhaps, one of possible ways 
to remedy this drawback is to gauge-rotate a Hamiltonian with non- 
polynomial eigenfunctions using a semi-classical or somehow modified 
semi-classical wavefunction as a gauge factor requiring an appearance 
of an algebraic form of the final operator. 
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